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We study the conductance of disordered wires with unitary symmetry focusing on the case 
in which m perfectly conducting channels are present due to the channel-number imbalance 
between two-propagating directions. Using the exact solution of the Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equation for transmission eigenvalues, we obtain the average and second 
moment of the conductance in the long-wire regime. For comparison, we employ the three- 
edge Chalker-Coddington model as the simplest example of channel-number-imbalanced sys- 
tems with m = 1, and obtain the average and second moment of the conductance by using 
a supersymmetry approach. We show that the result for the Chalker-Coddington model is 
identical to that obtained from the DMPK equation. 
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1. Introduction 

The statistical property of electron transport in a disordered quantum wire is indepen- 
dent of microscopic details of the system, and is mainly determined by the symmetries the 
system possesses. 1 In ordinary disordered systems, only time-reversal and spin-rotation sym- 
metries play a relevant role. According to the presence or absence of these two symmetries, 
disordered quantum wires are classified into either of the three standard universality classes 
(i.e., orthogonal, unitary and symplectic). The orthogonal class consists of systems having 
both time-reversal and spin-rotation symmetries, while the unitary class is characterized by 
the absence of time-reversal symmetry. The systems having time reversal-symmetry without 
spin-rotation symmetry belong to the symplectic class. It has been believed that although de- 
tails of transport properties differ from class to class, Anderson localization inevitably arises 
in all the three standard classes (i.e., the conductance decays exponentially with increasing 
sample length L and eventually vanishes in the limit of L — > oo). However, this widely ac- 
cepted understanding dose not always hold true. Recent studies on the symplectic universality 
class show that if the number of conducting channels is odd, one channel becomes perfectly 
conducting without backward scattering. 2-8 Due to the presence of this perfectly conducting 
channel, the dimensionless conductance behaves as g — * 1 with increasing L, and thereby An- 
derson localization disappears. In contrast to this, such an anomalous behavior does not arise 
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in the ordinary case with an even number of conducting channels. Thus, we must separate 
the symplectic class into two subclasses of symplectic-even and symplectic-odd according to 
whether the number of conducting channels is even or odd. 3 ' 7 

The realization of the symplectic-odd class revives attention to disordered wire systems 
with the channel-number imbalance between two-propagating directions. More than a decade 
ago, Barnes, Johnson and Kirczenow 9 pointed out that if the number of conducting channels in 
one propagating direction is by m greater than that in the opposite direction, the dimensionless 
conductances g for the majority direction and g' for the minority direction satisfy g = g' + m. 
Performing a numerical simulation, they obtained evidence that g — ► m and g' — > in the 
long-L limit. This indicates that m perfectly conducting channels are present only in the 
majority direction. Note that the channel-number imbalance leads to the breaking of time- 
reversal symmetry, so that their argument applies only to the unitary class. Recently, Hirose, 
Ohtuski and Slevin 10 have proposed the Chalker-Coddington model 11 with an odd-number of 
edge channels and studied its electron transport properties. The Chalker-Coddington model 
can be viewed as a stacking of alternating left-moving and right-moving chiral edge channels 
with random tunneling between adjacent ones. 12 It should be noted that in the odd-edge 
case, the channel number in one propagating direction is by one greater than that in the 
opposite direction. This situation is equivalent to that in the case of m = 1 considered in 
ref. 9. Performing a numerical simulation, they also confirmed that g — > 1 and g' — > in the 
long-L limit. 

One may think that the channel- number imbalance is rather unrealistic in actual systems. 
However, Wakabayashi et a/. 13 have recently shown that such a system can be realized in zigzag 
nanographene ribbons. 1 Inspired by this observation, the present author and Wakabayashi 15 
have formulated a random-matrix theory for the unitary universality class with the channel- 
number imbalance. Let us focus on a disordered wire system of length L having N + m left- 
moving channels and N right-moving channels. In this case, m left-moving channels become 
perfectly conducting and the dimensionless conductances g and g' for the left-moving and 
right-moving channels, respectively, satisfy g = g' + m. They have derived the Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equation, which describes the evolution of the probability 
distribution for transmission eigenvalues with increasing L, 16 ' 17 and analyzed the asymptotic 
behavior of g' = g — m in the long-L regime by using an approximate method proposed by 
Pichard. 18 They have shown that the localization length £, which characterizes the exponential 
decay of exp[(ln</)], depends on m as £ = 2Nl/(m + 1), where I is the mean free path for 
the left-moving channels. That is, £ decreases with increasing the number m of perfectly 
conducting channels. This means that the presence of perfectly conducting channels suppresses 
g'. Furthermore, they have also shown that the average and second moment of g' behave 
as (<?') ~ e~4jv for m = 0, (g'} ~ e~^~ for m > 1, (g' 2 } ~ e _( 4iv for < m < 3, 
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2 2(m-l)« 

and (g ) ~ e n for m > 3, where s = L/l. To examine these results in terms of 
an independent approach, the present author and Wakabayashi 19 have studied the averaged 
conductance (</) in the three-edge Chalker-Coddington model, i.e., the simplest nontrivial 
example of the channel-number-imbalanced unitary class with m = 1. They have found that 
(g') ~ e~^v, which is consistent with the DMPK result. Although this supports the validity 
of the DMPK approach, a more detailed comparison is highly desirable to deeply understand 
the nature of this peculiar universality class. 

In this paper, we study the average and second moment of the dimensionless conductance 
g' for the channel-number-imbalanced unitary class. First, we obtain the asymptotic forms 
of (g 1 ) and (g' 2 ) in the long-L regime for < m < 4 by using the existing exact solution 20 
of the DMPK equation. The exact solution, which is available only in the unitary class, 21 
enables us to obtain the full asymptotic forms including a pre-exponential factor. 22 Second, 
we employ the M-edge Chalker-Coddington model for the cases of M = 2 and M = 3, and 
obtain (g'} and (g' 2 } in the long-L regime by using a supersymmetry approach. 23 The M = 2 
case is the simplest example of the ordinary unitary class, while the M = 3 case is that of 
the channel-number-imbalanced unitary class. We obtain the full asymptotic forms of {g') 
and (g' 2 ) including a pre-exponential factor. We show that the resulting asymptotic forms are 
identical to those obtained from the DMPK equation. This strongly supports the validity of 
the DMPK approach. 

In the next section, we introduce the DMPK equation for the channel-number-imbalanced 
unitary class and introduce its exact solution. We obtain the average {g') and second moment 
(g' 2 ) of the dimensionless conductance g' in the long-L regime using the exact solution. In §3, 
we introduce the M-edge Chalker-Coddington model and obtain (g') and {g' 2 } for the cases 
of M = 2 and M = 3. Section 4 is devoted to summary. 

2. DMPK Approach 

We summarize the random-matrix theory for the channel-number-imbalanced unitary 
class. 15 We consider the case in which the number of left-moving channels is N + m, 
while that of right-moving channels is N. In this case, we can show that m transmis- 
sion eigenvalues in the left-moving channels become unity. 9, 10, 15 This indicates the pres- 
ence of perfectly conducting channels. If the set of the transmission eigenvalues for the 
right-moving channels is {Ti, T2, . . . , Tn}, that for the left-moving channels is expressed as 
{Ti, T2, . . . , Tjy, 1, . . . , 1}, where we have identified the N + 1 to N + mth channels as the per- 
fectly conducting ones. The dimensionless conductance g for the left-moving channels is given 
by g = S^i™ T a = m + J2a=i Ta, while that for the right-moving channels is g' = J2a=i ^a- 
It is easy to observe that g = g' + m. We consider the behavior of g' = g — m as a function 
of the normalized system length s = L/l, where / is the mean free path for the left-moving 
channels. It should be noted that the mean free path for the right-moving channels is not 
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equal to I due to the presence of perfectly conducting channels only in the left-moving chan- 
nels. Indeed, we find that I' = (N/(N + m))l. We define A a = (1 — T a )/T a and introduce 
the probability distribution P({A a }; s) for the transmission eigenvalues. The Fokker-Planck 
equation for P({X a }; s), which is usually called the DMPK equation, is expressed as 15 

with 

N N-l N 

j =ii a ™ x n n i^-^i 2 - (2) 

c=l 6=1 a=b+l 

As stressed in ref. 15, the factor Y\^=i m J represents the repulsion arising from the m- 
fold degenerate perfectly conducting eigenvalue. This reduces the non-perfectly conducting 
eigenvalues {Ti, T2, . . . , Tn}. It should be mentioned that the equivalent DMPK equation 
was proposed by Akuzawa and Wadati 20 in a rather formal context, but they did not study 
electron transport properties. 

The DMPK equation has been solved exactly for the ordinary case of m = 0. 21 ' 22 The 
exact solution for an arbitrary m has been obtained in ref. 20. In our notation, the probability 
distribution is given by 

P(M; s) = ^Bet{K(X a , X b ; s)} a , b=1 , 2 ,..., N (3) 



with 



Here, 



N 

K(\,X , ; 8 ) = ^Q j (\,8)h j (X , ,8). ( 4 ) 
j'=i 



(2j+m-l) 2 



Qj(X, s) = X m G j _ 1 (m + l,m + 1; -X)e sv — s } (5) 

r 

v,*)= / 

J 



hj(X,s)= dfcLj(fc)<4(fc)e~** a F , ,m + l;-A , (6) 



where Gj-\(m + l,m + 1; — A) = F(— j + l,m + j,m + 1; — A) is the Jacobi polynomial and 

fc 2 + (2l + m-l) 2 
-(2j+m-l) 2 + (2Z + m-l) 2 ' 



L(k)- TT fc 2 + (2/ + m-l) 2 



Cmlj "v^ r(m + i)|r(ijb)| • w 



Note that if (A, A'; s) satisfies 

dAiT(A, A; s) = N, (9) 



poo 



poo 



PC 

/ dA"K(A, A"; s)K(X", A'; s) = if (A, A'; s). (10) 
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Using these equations, we obtain the eigenvalue density 

poo 

i?i (An s) = N dA 2 • ■ ■ dAivP({A a }; s) 
Jo 

= K(X l ,X l ; S ), (11) 
and the two-point correlation function 

POO 

i? 2 (Ai, A 2 ; s) = N(N - 1) / dA 3 • • • dX N P({X a }; s) 

Jo 

= K{X 1 ,X 1 ;s)K(X 2 , A 2 ; s) - K(X 1 , A 2 ; s)K(X 2 , A l5 s). (12) 

We focus on the dimensionless conductance g' for the right-moving channels, in terms of which 
the dimensionless conductance g for the left-moving channels is given by g = m + g' . Using 
the eigenvalue density, we express the averaged dimensionless conductance as 

ti) = I dA— i?!(A;s). (13) 

Using the eigenvalue density and the two-point correlation function, we express the second 
moment as 

/•OO -I j-OO POO 1 

y 2 ) = j M (TTW R > {X:s) + l iX L dx \i + W + y) mx ' x '-' a) - (14) 

Alternatively, once (g 1 ) is given, (g' 2 ) is obtained by using the scaling relation 

N^l = -(g' 2 )-m(g'), (15) 

which can be derived from eq. (1) (see Appendix A). Using these equations, we obtain (g f ) and 
(g' 2 } in the asymptotic regime of s = L/l S> 4A^. It should be mentioned that the ordinary 
case of m = has been analyzed in ref. 22. 

We obtain the asymptotic form of (g'). We rewrite eq. (13) as 

f -oo i N 

Jo 1 + X J=1 

From eq. (5), we observe that the term with j = 1 dominates the others in the asymptotic 
regime of s S> AN, so that we can neglect the terms with j > 2. Equation (16) is then reduced 
to 

fc 2 + (m + l) 2 



poo \m roo fc 2 +( 

{9 '™ ) = J dX TTxJ dkL ^ c2 m( k > — ^ 



xFi^ , ,m + l;-Aj, (17) 

where Go(m + l,m + 1; —A) = 1 has been used. Here and hereafter, we explicitly show the 
number of perfectly conducting channels as a subscript. Before describing the evaluation of 
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eq. (17), we present the final results, 

3 

Wo) = ^ra N ——~e-^, (18) 



4 



AN 



2 1 

<s'i> = — a N^^e N , (19) 



7T2 ViiV 

W = ^Jve-£, (20) 
(5 3 ) = ^-a N e~^, (21) 



67T 

<04> = ^^"^ ( 22 ) 



where 



2 



Note that ai = tt/4 and liniTv^ooajv = 1- We observe that the exponential decay of (g' m ) 
becomes faster with increasing m. This should be attributed to the eigenvalue repulsion arising 
from the m-fold degenerate perfectly conducting eigenvalue. 

We briefly describe the derivation of eqs. (18)-(22). For the cases of m = and m = 1, we 
can exchange the order of the integrations over A and k in eq. (17) and find 

(g'J = / dkL^cK^e ^ s I m (k) (24) 

J 

with 

T „s f°° A m /m+l-ifc m + l + ik n A 
Adapting the method presented in ref. 24, we analytically obtain I m (k) as 

m = ^(fy (26) 



h{k) = t-tt^. (27) 



2tt 

/csimTff ) 

We next carry out the integration over k in eq. (24). The main contribution comes from the 



small-A: region of k < ^JAN/s, and therefore we approximate as L\(k) = Li(0), Cq(/c) = 

(7r/4)fc 2 , c 2 (/c) = (47r) _1 A; 2 , Io{k) = ir and ii(fc) = 4/fc 2 . After the fc-integration, we finally 

obtain eqs. (18) and (19). For m > 2, we employ a different approach, which is applicable to the 

cases of m > 1. We consider eq. (17) in the large-s limit. Expecting that the main contribution 

to the A-integration comes from the region of A ^ 1, we replace the hypergeometric function 

in eq. (17) by its asymptotic form 

m + 1 — ik m + l + ik A Tim + l)T(ik) . m+i-ifc , . 
, ,m + l;-A =— 9 A 2 +c.c. (28) 

2 2 J r ^m+l+ifc^2 v ) 
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in the large-A limit. We treat the case of m = 1 as an example. Note that the main contribution 
to the ^-integration comes from the region of k < \J 4iV/s, in which we can approximate as 
Li(k) = Li(0) and c\(k) = (47r) _1 fc 2 , and eq. (28) with m = 1 is reduced to 



2 + ifc ^ 2 _ _^ = 2A _ l S in (fc^) 



(29) 




^1) = ^ i-s e-^/ dA— -^exp I --( I. (30) 



After the fc-integration, we obtain 

MO) 

V V 4AT 

Changing the variable from A to x defined by A = (cosh(2x) — l)/2, we obtain 

r 1 lnA f iV/lnA\ 2 \ J" _* x 2 s , 01 . 

Substituting this into eq. (30) and using Li(0) = (4/7r)aAr, we again arrive at eq. (19). Adapt- 
ing the above method to the cases of m > 2, we obtain eqs. (20)-(22). 

We turn to the evaluation of the second moment in the asymptotic regime. The results 
are summarized as follows: 

(g'l) = ^a N -^e-^, (32) 



16 



AN 



{g'l) = -\a N ^=se n, (33) 

47T2 



AN 



3 1 

(5 2) = TT^aN—^e «, (34) 



128 



4Af 



<5l) = — "at— ^=e (35) 
3vr2 y 4]v 



{g'l) = T^ ajve " • ( 36 ) 

Again, we observe that the decay of (g'm) becomes faster with increasing m due to the 
eigenvalue repulsion from the m-fold degenerate perfectly conducting eigenvalue. 

We briefly present the derivation of eqs. (32)-(36). From eq. (15), the second moment is 
expressed in terms of the averaged conductance as 

(g' 2 J = -N^--m( 9 ' m ). (37) 

Substituting eqs. (18) and (19) into this, we straightforwardly obtain eqs. (32) and (33), 
respectively. However, the right-hand side of eq. (37) vanishes if we substitute eqs. (20)- (22). 
This means that the second moment for m > 2 is related to the next leading order correction 
to (g' m )- Instead of evaluating such a correction, we directly obtain using eq. (14). In 
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the asymptotic regime, eq. (14) is approximately reduced to 

P °° X" 1 f°° „ , s , ,, s fc 2 +(m + l) 2 



poo ym poo 

Wm)=J Q dA ( 1 + A ) 2 y o d/cLi(/c)c^(fc)e" 



4AT S 



— ifcm+l + i/c \ . . 

x F , ,m + 1; -Aj . (38) 

When m < 3, we can exchange the order of the A- and /c-integrations and find 

, 9 f 00 fc 2 +(m+l) 2 _ 

<</ 2 > = / dkL^cl^e — s I m {k) (39) 

J o 



with 



oc 



A m + 1 — i/c m + l + i/c 



/m(fc)= y dA (TTAF F i^^'^^' m+i; " A j- (4o) 

Adapting the method presented in ref. 24, we can analytically obtain I m (k) as 

m =d?m- (41) 
h{k)= i£wr <42> 

After carrying out the integration over k, we obtain eqs. (34) and (35). This approach cannot 
be applied when m > 4, so we adapt the method used to derive eqs. (20)-(22). It is applicable 
to the cases of m > 3. We replace the hypergeometric function in eq. (38) by its asymptotic 
form and integrate over k. For m = 4, we obtain 

,,2, ■i-^LiiO) 1 -a r ,. AI In A ( N (\n\\ 2 \ 



4Af 

After the A-integration, we finally arrive at eq. (36). 

3. Chalker-Coddington Model 

We consider an electron system consisting of M chiral edge channels as shown in Fig. 1. 
The propagating direction of the jth edge channel is left (right) if j is odd (even), and each 
edge is coupled with adjacent counter-propagating edge(s) by random tunneling. Note that 
this system is essentially equivalent to the Chalker-Coddington model. We assume that the 
electron wavefunction (f)j(x) for the jth edge channel with energy e obeys 12 ' 23 

(-iy(-i)d x <f>j(x) + tj(x)<l> j+1 (x) + ^(aO^-itx) = e^-(x), (44) 

where 4>o( x ) = (pM+iix) = has been assumed. Here, tj(x) and tj(x) represent the tunneling 
amplitude from the jth edge to j + 1th edge and that for the reverse process, respectively. We 
assume that {tj(x)} are random variables in the disordered region of < x < L, and vanish 
outside this region. We regard L as the length of our system, and the clean regions of x < 
and L < x play a role of the left and right electrodes, respectively. The random amplitudes 
in the disordered region are assumed to be Gaussian distributed with zero means. That is, 
(tj(x)) = 0, {tjixYj^x')) = and 

(t 3 (x)t*,(x')) = D5^d(x-x'), (45) 
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Fig. 1. Schematic sketch of the Chalker-Coddington model for (a) the two-edge case of M = 2 and 
(b) the three-edge case of M = 3. 



where (• • • ) represents the ensemble average. We hereafter set e = without loss of generality. 

The dimensionless conductance g for the left-moving channels is expressed as g = 
^ ij=odd 7ij, where %j is the transmission probability for an electron incoming from the 
jth channel in the right electrode and outgoing to the ith channel in the left electrode. In a 
manner similar to this, g' for the right-moving channels is expressed as g' = J2i j= C vcn ^ij- The 
dimensionless conductances satisfy g = g' in the ordinary case of an even M. In the odd-M 
case, which is first studied by Hirose, Ohtuski and Slevin, 10 the number of the left-moving 
channels is by one greater than that of the right-moving channels, so that one left-moving 
channel becomes perfectly conducting. In this case, we observe that g = 1 + g' . We consider 
the two-edge case of M = 2 and the three-edge case of M = 3 in the following. In the notation 
used in the previous section, the former and latter correspond to the cases of N = 1 with 
m = and N = 1 with m = 1, respectively. The three-edge case of M = 3 is the simplest 
nontrivial example of the disordered wires with the channel-number imbalance. Hereafter, we 
again explicitly show the number m of perfectly conducting channels as a subscript of g' . We 
focus on the average and second moment of g' m (m = 0, 1) as a function of L. Note that 
g' m = T22 in the two- and three-edge cases, in which only one right-moving channel exists. 

Mathur has shown for the M = 2 case that 23 ' 25 



I PL 

V 4 

in the asymptotic regime of DL » 1. This result is identical to eq. (18) in the case of N = 1 if 
we replace DL by s. This indicates that D" 1 should be identified with the mean free path I for 
the M = 2 case. In the M = 3 case, each left-moving edge channel is directly coupled to one 
right-moving edge channel as in the case of M = 2, so the mean free path for the left direction 
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is equal to I under the condition of eq. (45). However, because the right-moving channel 
interacts with the two left-moving channels, the corresponding mean free path becomes half 
of I. Thus, for comparison with the argument in §2, in which the system length is normalized 
by the mean free path for the left-moving channels, we hereafter adopt D~ l as the unit length 
scale common to both the cases of M = 2 and 3. This is simply achieved by the replacement 
DL -► s. 

We obtain the asymptotic forms of (g'^}, (g'i) and (g'\) by adapting the supersymmetry 
approach presented by Mathur. 23 The supersymmetry approach enables us to express the av- 
erage and second moment of the transmission probability in the from of a correlation function 
for a non-random interacting fermion-boson system on a one-dimensional chain. For the two- 
edge (three-edge) case, we need to consider the two-site (three-site) chain which we describe 
below. Let Ci a (cj CT ) be the fermion annihilation (creation) operator, and h%a (b\ a ) be the boson 
annihilation (creation) operator, where i = 1,2,3 represents the site number and a =|, j. The 
sites are one-to-one correspondence with the edge channels in the original model. The Hamil- 
tonian H for the two-edge case is identical to H\ 2 which represents the interaction between 
the 1st and 2nd sites. The Hamiltonian for the three-edge case is given by H = H\ 2 + H23, 
where H23 represents the interaction between the 2nd and 3rd sites. The explicit form of Hyi 
is given by #12 = Hf 2 + Hf 2 + Hf 2 B with 



Hf 2 = - ( c L c i^ - V 2 ) ( c L c 2a - I/ 2 ) + c^c^c^i + ci|Ci T c 2i c 2 t, (47) 
<r=T,l 

H * = E ( 6 i>- + V2) (bl a b 2a + 1/2) - b\^ n b\^ n - b n b n b 2i b 2h (48) 

<T=T,1 

H 12 = E (~ C L & l^ c L & 2a - C^bl^abla + Aa h \a c 2a h \a + CxM^Mci) , (49) 

where a in eq. (49) means f =[ and I =]. The simple replacement 1 — > 3 in the expression 
of H12 yields H 2 ^. In terms of the vacuum state |0), the average and second moment of the 
transmission probability are expressed as 23 

(T ij ) = (0\c ii c i ie- 6 ^J i ,\0), (50) 



(T?) = (OlciiCiT^e-^&U}, cU, |0) (51) 



with Dl = DL. Here, c^Cjf, bifi^ and their Hermitian conjugates play a role of the current 
vertex, where j and J, correspond to the retarded and advanced sectors, respectively. We see 
that H is non-Hermitian because Hf 2 B = —Hj 2 and = —H 2 B , but it has only real 

eigenvalues. Equation (51) has not been presented in ref. 23, but we can easily derive it by 
extending Mathur's argument. 

Our task is now reduced to evaluating the correlation functions for the non-random in- 
teracting fermion-boson system. Note that Hf 2 B and H 2 B annihilate the two-fermion state 
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c^cljO), while iJf 2 c 2-f c 2j.l^) = ^23 C 2T C 2|I^) = (V^c^c^JO). Consequently, the expressions 
for (g' m ) = (T22) and (g'm) = (^22) are simplified to 

<«/ m ) = (0|e-^|0), (52) 

W 2 J = <0|6 2 |&2Te-^4AlO> (53) 
with Hff = Hf 2 + 1/2 and Hf = Hf 2 + H\\ + 1. Equations (52) and (53) indicate that we can 
evaluate (g' m ) and (g' 2 m ) in the asymptotic regime of L D^ 1 if the low-energy eigenstates 
of the boson Hamiltonian H% (m = 1,2) are obtained. That is, we need not consider the 
fermion degrees of freedom in the following argument. This enables us to treat our task in 
an analytical manner. It should be emphasized that this special simplification arises only in 
the case that we treat the right-moving conductance g' with M = 2 or 3. In other words, an 
analytical treatment seems to be difficult for M > 4. 

To evaluate (g'l), we consider eigenstates |^) of Hq. We need to obtain the eigenstates 
having an overlap with an d therefore we focus on the Hilbert space spanned by 

i^ n ) = ^Vt)i ( b ^T ovy n+1 1°) ™ 

with n = 0, 1, 2, . . . , 00. Eigenstates of Hq is expressed as |^) = J2T=o c n\4> n }- From the 
eigenvalue equation Hq\^} = E\^}, we obtain the recurrence relation 

(2n 2 + 4n + 2 - E)c n - n(n + l)c n _i - (n + l)(n + 2)c n+1 = 0. (55) 

On the basis of eq. (55), we obtain the continuous eigenstates \k) which satisfy H®\k) = 
\{1 + k 2 )\k) and (0| 624, I ^) = 1 w ^h k > 0. The orthogonality condition is 

The derivation of the above result is briefly described in Appendix B. From eq. (56), we 
observe that the identity operator X is expressed as 



L 



00 7r/csinh(^) , w , 

dk ±*±\k){k\. 57 

8cosh(^) 1 n 1 ^ J 



We now evaluate (g'o)- Inserting eq. (57) into eq. (53), we obtain 

The main contribution to the /c-integration comes from the region of k <C 1, in which we can 



approximate as 



8 cosh 16 

After the /c-integration, we obtain 



7rfcsmh(f)_^ ^ 



7T2 1 



DL 



DL 
4 
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This result is identical to eq. (32) in the case oi N = 1 under the replacement DL — > s. 

To evaluate (g'i) in the asymptotic regime, we consider low- lying eigenstates \^>) of Hf 
having an overlap with |0). 19 Therefore, we focus on the Hilbert space spanned by 

= „.(„+„,))„,! ( W ( b lAiT +m (*A) m m m 

with n, m = 0, 1, 2, . . . , oo. Eigenstates of is expressed as = ^^° m=0 c ra ' m |V> n ' m ). From 
the eigenvalue equation Hf |*) = E 1 )^), we obtain the recurrence relation 

{1 + 2(n + m + l)(n + m + 1/2) - £} c n ' m - n(n + m)c"- 1>m - (n + r^mc™'™- 1 

- (n + l)(n + m + l)c n+1 ' m - (n + m + l)(m + l)c n ' m+1 = 0. (62) 

We restrict our attention to the lowest energy branch of the excitation spectrum, which 
determines the asymptotic behavior of (g\). We assume for this branch that c n ' m depends on 
only n + m (i.e., c n+m '° = • • • = c n,m = • • • = c o,n+m^ which h as been confirmed in ref. 19 by 
a numerical diagonalization of eq. (62). Under this assumption, we rewrite c n,m as c n,m — »■ q 
with I = n + m. Equation (62) is then reduced to 

{1 + 2(1 + 1)(Z + 1/2) - £} q - / 2 q_i - (/ + 1)(Z + 2) Q+1 = 0. (63) 

On the basis of eq. (63), we obtain the continuous eigenstates \k) which satisfy Hf\k) = 
(1 + k 2 )\k) and (0\k) = 1 with k > 0. The orthogonality condition is 

<*!*>- i£^S)' (*-*>■ < 64 > 

The derivation of the above result is briefly described in Appendix C. From eq. (64), we 
observe that the identity operator 1 in the restricted Hilbert space is expressed as 

27rfeWi(7rfe) 
Sinn (7rfe) 

Inserting eq. (65) into eq. (52), we obtain 

W-f ^fflf K°WI 2 ^" +fcg) - <«» 

JO Smh (7TK) 

Carrying out the /c-integration, we obtain 

2^/7T / DL 
V 4 

This result is identical to eq. (19) in the case oi N = 1 under the replacement DL — > s. 

Finally, we evaluate (</i) in the asymptotic regime. To do so, we consider low- lying eigen- 
states |3>) of Hf having an overlap with frJ^JO). Therefore, we focus on the Hilbert space 
spanned by 

n+m+l 



f 

JO 



n\(n + m + l)!m! 



( &t i T fet u) n ( 6 vy n m (4 T & y m io) (es) 



12/20 



J. Phys. Soc. Jpn. Full Paper 

with n, m = 0, 1, 2, . . . , oo. Eigenstates of Hf are expressed as |$) = Y^nm=o d n,m \(j) n,m ) . From 
the eigenvalue equation Hf\<&) = E\$>), we obtain the recurrence relation 

{1 + 2(n + m + l)(n + m + 3/2) - E} d n ' m - n(n + m + l)^" 1 '" 1 - (n + m + ljmff-^ 1 

- (n + l)(n + m + 2)d n+1 ' m - (n + m + 2)(m + l)d n ' m+1 = 0. (69) 

Again, our attention is restricted to the lowest energy branch of the excitation spectrum. 
We assume for this branch that d n ' m depends on only n + m (i.e., d n+m '° = • • • = d n,m = 
• • • = d 0,n+m ), which can be confirmed by a numerical diagonalization of eq. (69). Under this 
assumption, we rewrite d n,m as d n,m — > with I = n + m. Equation (69) is then reduced to 

{1 + 2(1 + 1)(Z + 3/2) - £} ^ - /(/ + l)d,_! - (/ + 2) 2 d i+1 = 0. (70) 

On the basis of eq. (70), we obtain the continuous eigenstates \k) which satisfy Hf \k) = 
(1 + k 2 )\k) and (0|62|&2tl^) = — k 2 with k > 0. The orthogonality condition is 

<™ = sSy' <*-*>■ («) 

The derivation of the above result is briefly described in Appendix C. From eq. (71), we 
observe that the identity operator J in the restricted Hilbert space is expressed as 

sinn (irk) 

Inserting eq. (72) into eq. (53), we obtain 

/ ' 2 \ f°° 27r cosh(7rfc) 2 _D L (l+fc 2 ) / VQ A 

(9i)= dfc |(0|6 2 |&2Tl fc )l e • ( 73 ) 

jo sinn (irk) 

Carrying out the /c-integration, we obtain 



poo 

=1 ' 



4 

This result is identical to eq. (33) in the case of N = 1 under the replacement DL — ► s. 
4. Summary 

We have studied electron transport properties in disordered unitary wires of length L 
in the presence of the channel-number imbalance between two propagating directions. Our 
attention is focused on the case in which the number of left-moving channels is by m greater 
than that of the right-moving ones. In this case, m left-moving channels become perfectly 
conducting and the dimensionless conductances g and g' for the left-moving and right-moving 
channels, respectively, satisfy g = g' + m. First, we have obtained the average (g') and second 
moment (g' 2 ) of g' = g — m in the long-L regime by using the exact solution of the DMPK 
equation. Both (</) and (g' 2 } decay exponentially as a function of L. It is shown that their 
exponential decay becomes faster with increasing m. This behavior can be understood from the 
fact that the eigenvalue repulsion arising from the perfectly conducting eigenvalue is enhanced 
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with increasing m. Second, we have employed the M-edge Chalker-Coddington model for the 
cases of M = 2 and M = 3, and obtained (g 1 ) and {g' 2 ) in the long-L regime by using the 
supersymmetry approach. The case of M = 3 corresponds to the simplest nontrivial example 
of the channel-number-imbalanced unitary class. We have shown that the resulting asymptotic 
forms of (g 1 ) and (g' 2 ) are identical to those obtained from the DMPK equation including the 
pre-exponential factor. 

Appendix A: Derivation of the Scaling Relation 

The average of a function F({A a }) is defined as 

/•oo 

{F)= d\ 1 ---X N F({X a })P({X a },s). (A-l) 
J o 

Using the DMPK equation given in eq. (1), we find that 26 



\a=l 

N 



( r) 2 F r)F 

g| Aa(1 + ,j_ + (1 + 2Aa) _ 



(b#o) 

Replacing F by g' = X)£Li(l + -^a)" 1 i n t ne above equation, we obtain eq. (15) after straight- 
forward calculations. 

Appendix B: Energy Spectrum of 

As shown in the text, eigenstates of Hq satisfying (0|&2J.&2tI*) / are expressed as 

oo 

\*{E)) = Y J C n {E)\r), (B-l) 

n=0 

where c n satisfies eq. (55). We obtain the asymptotic form of c n in the large-n limit adapting 

the manipulation given by Mathur. 23 We introduce the generating function defined by f(x) = 

Y^=o c n xn - Using eq. (55), we can show that f(x) obeys 

d 2 f df 
x(x - l) 2 -4 + 2{2x -l)(x- + {2x - 2 + E)f = 0. (B-2) 

dx z dx 

It is convenient to rewrite f(x) as f(x) = (1 — x)^g(x). If we set 

V = ~\ + K (B-3) 

with k = iy/AE - 1/2 for E > 1/4 and k = y/l - 4E/2 for E 1 < 1/4, then is expressed 
in terms of the hyper geometric function as g(x) = F(2 + n, 1 + /j,, 2; x). The coefficient c n is 
expressed as 

c n = ±-J c dxJ ri (l-xr 9 (x), (B-4) 
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where C denotes a small contour encircling the origin in the anticlockwise direction and cq = 1 
has been assumed. We employ the integral representation of the hypergeometric function 

F(a, b, c; x) = TQ> ^_ fe) J™ dt(t - - l)*"*" 1 , (B-5) 

which is justified when Re[c] > Rc[b] > 0. Substituting this into eq. (B-4) and exchanging the 
order of the integrations over x and t, we obtain 

1 r°° /■ rt.r I 

Cn -r(i + /,)r(i-/,)A 

We here draw the branch cut between x = 1 and x = t, and take the phase of (t — x) and 
(1 — x) being zero when x lies on the real axis to the left of 1. Deforming the contour C as 
in ref. 23, we can show that C is replaced by the contour starting from x = 1 to x = t above 
the branch cut and after encircling the point x = t in the clockwise direction, coming back to 
x = 1 below the branch cut. Exchanging the order of the integrations, we obtain 

- " - r ( i + ,)r ( i-M) f dl ^/ c , £ ( " (( " W ~ ^ (B - 7) 

where we have drawn the branch cut between t = x and t = oo on the real axis and C x denotes 
the contour coming from t = oo to t = x below the branch cut and after encircling the point 
t = x in the clockwise direction, going back to t = oo above the branch cut. Rescaling t as 
u = t/x, we obtain 

Cn = -i- / dx ^— 4 a(x) (B-8) 

r(i + /x)r(i - n) J 1 x n + 2 +» K ' v ' 



with 



a(x)= [ |?V( u -ip»(l- u )-2-M (B-9) 



where Ci is identical to C x with a; — > 1. Using an analytic continuation, we can show that 



a(x) = —F ( 1 + fi,2 + fx,2;— ) . (B-10) 



Then, we obtain 



oo 



where /x = — 1/2 + k has been used. The change of the variable from x to s = lux results in 
c n = 1 f°° ds(e° - l)-^e-( n+ l^>F (\ + «, § + K , 2; e" s ) . (B-12) 

r(i + K)r(i - k) Jo \2 2 y 

The above equation indicates that the behavior of the integrand for s « is important in 
considering the large-n limit. Therefore, we employ the approximations (e s — 1)~2+ K m s _ 2+ K 
and 

Fin h r- P - s \ ~ r ( C ) r ( C ~ Q - fe ) , r(c)r(Q + 6 - C) s , c _ a _fe , , 

F(a ' 6 ' c ' e ^r^-ajrcc-fr) 4 " r(a)r(6) (1 " e } (B " 13) 
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which are applicable when e~ s ~ 1. Carrying out the s-integration, we obtain 

1 ( r(2K)r(l- K ) e Kln " T(-2k)T$ + k) e -^n \ 

If E 1 > 1/4, we can set k = ik/2 with = y/4E — 1. In this case, c n in the large-n limit is 
expressed as 

c n (£) = ^cosQlnn + C fe ) (B-15) 

with 

4 cos h 

"fc = 2 , (B-16) 

7ry fcsinh(7r/c) 

/ r(ifc)r(i-i§) \ 

Cfe = arg / V 2 „ 2 \ . B-17 

We show that the orthonormalization of \^f(E)) is possible when E > 1/4. Note that 
{^(E')\ i lf(E)) = Y^=o c n(E')c n {E). Using eq. (55), we can express the partial sum 
En= Cn(E')c n (E) as 

J2c n (E')c n (E) = ( M ^K^ + 2 ) ( CM+1 (£') CM (£) _ cm(E')c m+ i(E)) . (B-18) 

n=0 

Substituting eq. (B-15) into eq. (B-18) and then taking the limit of M — > oo, we obtain 

oo 

c n (E')cn(E) = ira 2 k 5 (k - k') . (B-19) 

n=0 

This indicates that \^>{E)) can be orthonormalized when E > 1/4. In contrast, when E < 1/4, 
the partial sum J2n=o c n{E')c n {E) does not converges in the large-M limit even if E ^ E'. 
Hence, the orthonormalization is impossible. We conclude that \^(E)) exits only when E > 
1/4. It is convenient to use k instead of E. We rewrite the eigenstate as \^(E)} — > \k) with 
E = \{1 + k 2 ). Equation (B-19) is then rewritten as 

8 cosh (—) 

w-;*^' (*-*>■ (B ' 20) 

Finally, we note that cq = 1 means (0|&2J.&2tI&) = 1- 

Appendix C: Low Energy Spectrum of 

We first consider the low-lying eigenstates of Hf satisfying (0|\E') ^ 0. Our attention is 
restricted to the lowest energy branch in which the eigenstates are expressed as 

oo 

mm) = E c n+m (E)\rn, (c-i) 

n,m=0 

where q with I = n + m satisfies eq. (63). We obtain the asymptotic form of q in the large-/ 
limit. We introduce the generating function defined by f(x) = Yli^o c l xl - Using eq. (63), we 
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can show that f(x) obeys 

x{x - 1) 2 + (3x - 2)(x - 1)^ + (x - 2 + E)f = 0. (C-2) 



It is convenient to rewrite f{x) as f(x) = (1 — x) fl g(x). If we set = WE — 1 for > 1 and 
= v 7 ! — -E 1 for E 1 < 1, then g(x) is expressed in terms of the hyper geometric function as 
g{x) = F(l + /j,, 1 + //, 2;x). The coefficient q is expressed as 

q = 2_jTdx 5 i r (l-x)^(x) > (C3) 

where C denotes a small contour encircling the origin and Co = 1 has been assumed. Using 
eq. (B-5) and adapting the procedure described in Appendix B, we obtain 

sin(7r / u) f°° (x — 1) M / 1\ t . 

After changing the variable from x to s = Ins, we replace the hypergeometric function by the 
approximate expression given in eq. (B-13). Carrying out the s-integration, we obtain 

_ sin(^) / r(l - n)T(2fi) e ^ ln (' +1 ) T(l + / u)r(-2 / u) e ^ ln (' +1 ) \ 



When E > 1, we write fj, = ik with k = y/E — 1. In this case, q in the large-/ limit is 
expressed as 

ci(E) = ^cos(k\n{l + l)+r ]k ) (C-6) 

with 

sinh(Trfc) 

Pk - - , =, (C-7) 
7r/cy cosh(7rKj 

/ r(i-ifc)r(2ife) \ 

^ = arg U(l + ifc)r(iA:) J" (C - 8) 
We show that the orthonormalization of \ty(E)) is possible when E > 1. Note that 

oo 

= c n ' m {E')c n ' m {E) 

n,m=0 

oo 



^ + 1)q(£')q(£)- (C-9) 



z=o 

Using eq. (63), we can express the partial sum YlbLott + ^) c i{E')ci(E) as 

£(J + l) Cl (E') Cl {E) = ( M + V 2 jM + 2) ( Cm+i{e , )cm{e) _ Cm{E > )cm+i{E) ) . (C-10) 

1=0 

Substituting eq. (C-6) into eq. (C-10) and then taking the limit of M — > oo, we obtain 

£(J + IME'ME) = ^5(k- k>) . (C-ll) 
i=o 

This indicates that \^(E)) can be orthonormalized when E > 1. In contrast, when < 1, we 
can show that the partial sum diverges in the limit of M — > oo even if E E' , and thereby 
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the orthonormalization is impossible. We conclude that \^f(E)) exists only when E > 1. It 
is convenient to rewrite the eigenstate as \ty(E)) — > \k) with E = 1 + k 2 . Equation (C-ll) is 
rewritten as 

< v '*>- a^a&) ' (*-*>• (c ' 12) 

We obtain (0|fe) = 1 from Co = 1. 

We next consider the low-lying eigenstates of Hf satisfying (O^^tI^) / 0- We restrict 
our attention to the lowest energy branch in which the eigenstates are expressed as 

oo 

ME)) = £ d n+m (EW m ), (C-13) 

n,m=0 

where d\ with / = n + m satisfies eq. (70). We introduce the generating function defined by 
f( x ) = Yli^o^ix 1 . Using eq. (70), we can show that f(x) obeys 

x 2 {x - 1) 2 ^4 + x(x - l)(4x - 3)^- + {2x 2 + (E - 4)x + 1}/ = 0. (C-14) 

dx z dx 

It is convenient to rewrite f(x) as f(x) = (1 — x^x -1 *?^). If we set y = \\JE — 1 for E > 1 
and /i = \/l — E for E 1 < 1, then g(x) is expressed in terms of the hypergeometric function 
as g(x) = F(l + /j,,fi, l;x). The coefficient q is expressed as 

d i = T~- I dx-ij-a-x)^, (C-15) 

Z7T1 Jc* x x 

where C denotes a small contour encircling the origin. It should be noted that eq. (C-15) 
indicates do = /i 2 . For g(x), we cannot directly employ the integral representation of the 
hypergeometric function given in eq. (B-5) when [i is pure imaginary, so we modify g(x) as 

1 - u 2 

g(x) =F(l + fi,l + ft, 2; x) j^xF(2 + n, 1 + n, 3; x). (C-16) 

Now, we can apply eq. (B-5) to each term in the right-hand side of the above equation. 
Adapting the procedure described in Appendix B, we obtain 



, sin(^) [°° , ( ix-lY , 1 



+ 1 1 + x)^F+5fF { F (- 1 + 2 + M, 1; i) - F (-1 + M , 1 + c, 1; i) } j ■ (C-17) 

After changing the variable from x to s = lnx, we replace the hypergeometric functions by 
the approximate expression given in eq. (B-13). Carrying out the s-integration, we obtain 
, _ sin(7r/x) / r(l + 2/i)r(l - y) e ^ ln (m) r(l - 2y)F(l + y) e "f M'+i) \ 

di - ~^r~ ^ faT^ rn fxr^P J ' (C ' 18) 

where several higher order terms with respect to y are neglected. 

When E > 1, we write y = ik with = \/E — 1. In this case, d\ in the large-£ limit is 
expressed as 

dl(E) = T ^jsm(kln(l + l) + e k ) (C-19) 



18/20 



J. Phys. Soc. Jpn. Full Paper 

with 

sinh(^fc) 

7fc 1 =, n » ( c - 2 °) 

7ry cosh(7rfcj 



^ = ^ ( y i;: : ^ 1 ) • (c-2i> 



T(i + 2ifc)r(i -ifc) 
r(i + ife) 

We show that the orthonormalization of | is possible when E > 1. Note that 
= £~ (/ + l)di{E')di{E). Using eq. (70), we can show 

+ miE'^E) = ( M +^K^ + 2 ) 2 (d M+1 ( E ')d M {E) - d M (E')d M+1 (E)) . (C-22) 

2=0 



Substituting eq. (C-19) into eq. (C-22) and then taking the limit of M — ► oo, we obtain 

r 2 
2" 



£(J + = ^ {k ~ k') . (C-23) 

1=0 

This indicates that |<1?(-E)) can be orthonormalized when E > 1. In contrast, when E 1 < 1, we 
can show that the partial sum diverges in the limit of M — > 00 even if E ^ E' , and thereby 
the orthonormalization is impossible. We conclude that \$>(E)) exits only when E > 1. It is 
convenient to rewrite the eigenstate as \3>{E)) — > with £7 = 1 + A; 2 . Equation (C-23) is 
rewritten as 

We obtain (0|&2J.&2tI&) = —k 2 from do = fi 2 . 
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